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Interesting websites: 

http://www.fachgruppe-computeralgebra.de/cms/tiki-index.php?page=Rundbrief 

You can download the „Rundbriefe“ (mostly in German) starting in December 1987. All docu- 
ments are in pdf- format. This provides a very interesting history of development and spreading of 
CAS. The Rundbriefe contain plenty of information about the various computer algebra systems. 

The next sites are recommended by Prof, de Villiers. Many tha nks for his wonderful and inspiring 
Newsletter. 

I'll give a selection of his recommendations. You can find much more starting with his personal 
website: 


http://mysite.mweb.co.za/residents/profmd/homepage4.html 


You can download 2 volumes (together nearly 600 pages!) containing the complete proceedings 
of the ICMI conference (Taiwan May 2009): 

ICMI Study 19: Proof and Proving in Mathematics Education 

from: http://ocs.library.utoronto.ca/index.php/icmi/8 


Goto the website of the Mathematical Association of America and download among others an 
article on Teaching and Learning Differential Equations by Chris Rasmussen and Karen White- 
head: 


http://www.maa.org/t_and_l/index.html 

http://www.maa.org/t_and_l/sampler/rs_7.html 


The next website recommended offers plenty of geometric models and how to create them with 
paper and scissors: 

http://www.korthalsaltes.com/ 

http://www.korthalsaltes.com/three_pyramides_in_a_cube.htm 

Download contributions published in the Journal of Mathematics Education (most of them sub- 
mitted by Chinese researchers: 


http://educationforatoz.com/journalandmagazines.html 
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Dear DUS Members, 


I apologize for being late with DNL#74. At the end of June and begin of July I 
attended some interesting conferences (ACA09, CAD&ME and others) which needed 
some extra preparations. Then some clarifications were necessary to finalize the 
DEQME-contribution of this DNL. And as I wrote in my short info our webmaster 
Walter Wegscheider enjoys his very well deserved holidays. When he will be back 
from his cruise on the Mediterranean Sea he will upload this issue. I hope that it 
was worth to wait one month: 


ACA09 (Applications of Computer Algebra 2009) was an excellent organized con- 
ference which was held in Montreal. Many thanks to Kathleen Pineau, Michel Beaudin 
and Sides Picard who made the conference a full success. There were among others 
a rich Educational Session and a special DERIVE session. You can find the abstracts 
of the DERIVE-related contributions. DNL#75 will present the full papers. 

By the way, don't forget to mark July 2010 in your agenda. 

TIME 2010 will be held in Malaga, Spain from 6-10 July. We 
will have four extraordinary keynote speakers: 

Barbel Barzel, Germany 
Michel Beaudin, Canada 
Colette Laborde, France 


Eugenio Roanes, Spain 

The Conference website is www.time2010.uma.es. 


This DNL contains a variety of articles. You can find another of R. Schroder's pro- 
jects for the classroom, an interesting paper on the use of the TI-NspireCAS, one 
of my exercising programs (originally written for my granddaughter), a great trian- 
gle problem where you can use various CASs, the first part of a review of Nils 
Hahnfeld's Differential Equations' package for the handheld and last but not least 
informs Johann Wiesenbauer in his Titibits 37 (!!!) about other methods for factor- 
izing integers. 


So I do hope that everybody will find some useful, interesting or just delighting 
pages. I wish you all a pretty summer and I am looking forward to meeting you again 
in fall. 



Download all DAX-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be wannly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS -TI Newsletter will be. 


Next issue: September 2009 

Deadline 15 August 2009 

Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Luke-Rosendahi, GER 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Simulating a Graphing Calculator in DERIVE , J. Bohm 
Henon, Mira, Gumowski & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - interesting features, T. Roller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the TI, M. R. Phillips 

Advanced Regression Routines for the TIs, M. R. Phillips 

Where oh Where is IT? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ 

Snail-shells, Piotr Trebisz, GER 

A Conics-Explorer, J. Bohm, AUT 

Coding Theory for the Classroom?, J. Bohm, AUT 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Runge-Kutta Unvealed, J. Bohm, AUT 

The Horror Octahedron, W. Alvermann, GER 

and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS -TI User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-06604070480 

e-mail: nojo.boehm@pgv.at 
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Derive- and CAS-TI-User Forum 



Tom Barbara, Oregon, USA 

Hi Josef, 

Tha nks for adding me to the group. 

I am a long time user of Derive (from Version 1). However, I did not use it that much in recent years, 
although I did have the first windows version of the software. I recently purchased the latest version 
that is available. Besides the general algebraic tools I used derive’s matrix capabilities for problems in 
finite dimensional quantum theory as applied to Nuclear Magnetic Resonance. Back then 1 lobbied the 
Derive folks to include a function for taking the Kronecker Product of matrices, 

And they were happy to comply. However, I now find that this feature is no longer supplied with the 
software. I believe this came as an additional operator that one could load from a library and I guess it 
feel off the list of those that were included. 

Perhaps some of the users recall this function and still have a copy of the library it came in? Or per- 
haps someone has invented their own operator or knows how to use the operators that are supplied to 
write a Kronecker Product function. Back then, it did not appear to be possible with the current tools 
and that is why I suggested to Derive that they supply one. 

Thanks, 

Tom B 


Dear Tom. 

I didn't find the Kronecker product among my many DERIVE files. 

And I must admit that I didn't know about this special matrix product, but I could find its defi- 
nition. 

So I took the challenge trying to produce the respective function (see the attached file). 

I hope that this works properly. Then I'll include a respective note in the next DERIVE News- 
letter. 

kron_prod(a , b) := APPEND(VECTOR(APPEND_COLUMNS((VECTOR(VECTOR(a -b, k, DIM(a’)), 

#1: i , k 

i, DIM(a))) V ], DIM(a)Y| 


r i 2 

knon_prod 3 4 


5 6 7 

8 9 10 

11 12 13 
12 3 


kron_pnod [[10, 20, 30]], 4 

^ L 5 


5 

6 

7 

10 

12 

14 

8 

9 

10 

16 

IS 

20 

11 

12 

13 

22 

24 

26 

1 

2 

3 

2 

4 

6 

15 

18 

21 

20 

24 

28 

24 

27 

30 

32 

36 

40 

33 

36 

39 

44 

48 

52 

3 

6 

9 

4 

8 

12 

25 

30 

35 

30 

36 

42 

40 

45 

50 

48 

54 

60 

55 

60 

65 

66 

72 

78 


30 

60 

90 

40 

80 

120 

50 

100 

150 


5 10 15 6 12 18 
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Tom Barbara, Oregon, USA 

Hi Josef, 

Well your routine looks correct. Studying your example will be a good lesson in programming for me. 
Thanks so much. 


Fred J. Tydeman 

I notice that in a function definition, one can have multiple statements on one line, with a comma as the separa- 
tor between statements. 

However, outside of a function, when I try either a comma or a semi-colon to separate multiple statements on 
one line, they are considered errors. 

So, how does one do multiple statements on one line? 


Danny Ross Lunsford 

Use a PROG statement to group them, e.g. 

PROG(x := 10, y := SIN(x), DISPLAY ("Hello world")) 

This evaluates to whatever is returned by the last statement in the group. 

-drl 

Fred J. Tydeman 

I created a file a.mth with the contents of: 

PROG (x: =3, y : =4 ) 

[a : =1 , b : =5 ] 
x+a 

I did Simplify, Basic the x+a line and got 

x + 1 

Johann Wiesenbauer 

Hi Fred, 

Maybe I'm wrong, but it appears to me that you didn't understand the crucial difference between [x:=3, 
y:=4] and prog (x : =3 , y : =4 ) . In the first case, x and y will have the corresponding values after 
INPUTTING that line, in the second case x and y will have those values only after SIMPLIFYING that line. 
After all, prog (x : =3 , y : =4 ) is a program, as its name says, and programs will never be started automati- 
cally in input mode. Hope this helps, 

Johann 

Germain Labonte [grmlabonte@yahoo.ca] 

Hello, 

I've come accross your contact information from the Derive User Group web site. 

I've had some interest in muLISP (on which Derive was programmed). From what I can gather, the 
latest version of muLISP which was sold up to the first half of 2005 (along with Derive for DOS) in- 
cluded both the 16-bit and 32-bit kernels. 

Is this correct? If so, what would be the version number. I'm hoping that one day a copy of this ver- 
sion of muLISP will come on the "used-software" market (e.g. eBay) and I would like to know what 
questions to ask the seller. 

Sincere regards, 

Germain Labonte, Mississauga, Ontario, Canada 

Is there anybody, who can help? 


More User Forum on page 41 
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The Common Measure 

Roland Schroder, Celle, Germany 

Two wooden sticks of lengths a = 27 cm and b = 21 cm shall be sawed up - without using 
any measuring instrument - to as many pieces of equal length as possible. We can do this 
by putting the sticks flush together and separate the overhang 
c = a - b from the longer stick. Now we choose the two shortest remaining sticks and repeat 
the procedure. This describes a recursive algorithm which ends if we cannot choose the two 
shorter sticks because they all are of the same length (which is 3 in our case). 

a = 27 i c = a - b = 6 



b = 21 

d = b - c = 15 


e = d - c = 9 



f = e - c = 3 





, g = c - f = 3 


Common Measure = 3 

We see that 3 is the GCD (Greatest Common Divisor) of 27 and 21 and that this recursive 
procedure is nothing else than the Euclidean Algorithm. I'd recommend giving the pupils 
pairs of straws or paper stripes and letting them perform this activity. They should find out 
that this is a way to find the GCD of two numbers. 

Maybe that some of you prefer using the computer instead of a saw, then you can write in 
DERIVE: 


'27 21 ' 

21 6 
6 15 


ITERATED [b, |a 


b|], [a, b], [27, 21], 7) 


15 9 

9 6 


6 5 


3 3 


3 0 J 


Calculation ends when comparing the two shortest sticks the overhang o = 0. The common 
measure is the minimum positive number appearing in the procedure. 

Of course, we don't know that only 7 iteration steps are necessary. What happens if we do 
not enter the number of steps? 
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We transfer this recursion to DERIVE: 

ITERATES([b, a - 2-b], [a, b], [1, J2 - 1], 8) 


and we receive 


1 

72-1 

72-1 

3 - 2-72 

UJ 

1 

ro 

5-72 - 7 

5-72 - 7 

17 - 12-72 

17 - 12-72 

29-72 - 41 

29-72 - 41 

99 - 70-72 

99 - 70-72 

169-72 - 239 

169-72 - 239 

577 - 408 - 72 

577 - 408-72 

985-72 - 1393 


No expression is appearing twice. Approximating the matrix shows that all expressions are 
positive. This confirms - at least for the first 9 steps - that each expression can be sub- 
tracted twice from its predecessor. 

We would like to find the explicit representation of this sequence. For this purpose we use a 
method which can be taken as the standard one. We conjecture that it might be a geometric 
sequence (aq n ) ne N- So we substitute for the respective elements of the sequence in the re- 
cursion formula: 

aq n+1 =a-q n - 1 -2a-q n . 

After division by aq n " 1 we get q 2 = 1 - 2q which is a quadratic equation with solutions 
q 1/2 = ±V2 - 1 . As all elements of the sequence are positive - and they shall remain positive 
- we accept only the positive solution. Using the positive quotient we can define a geometric 
sequence: 


VECTOR (Ln, 


0/2 - l) n J, 


0 , 8 ') = 


0 1 


1 

72 

- 1 

2 

3 - 

2-72 

3 

5-72 

- 7 

4 

17 - 

12-72 

5 

29-72 

- 41 

6 

99 - 

70-72 

7 

169 ■ 72 

- 239 

8 

577 - 

408-72 . 


The right column of this table corresponds with the left column of the result of the recursion 
from above. 





p 
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So it seems to be clear that alternative removement starting with two sticks of lengths 1 and 
a/2 is leading to a geometric sequence of stick lengths. Because of q < 1 it will cnverge to- 
wards zero but never reach zero. The procedure for finding a common measure will never 
come to an end. There good reasons to suspect, that 1 and a/2 don't have a common meas- 
ure. We can name 1 and a/2 as incommensureable. 

The procedure can be used to create approximating representations as fractions together 
with a sound error estimation. DERIVE gives: 

( a /2 - 1 ) 20 = 0,0000000221 

or: 

22619537 - 15994428-V2 « 0,0000000221 . 

Solving for a/2 results in the difference of a fraction (approximation for a/2 represented as a 
fraction) and a number less 10 “ 14 (error estimation). 


The foilwing screen shots show the procedure performed on the Voyage 200 (in Sequence 
Mode): 


ITi-THOT T F]: 

IZoonlEdit 

FH T FS^Y FbT Y F7 f 'I 

^ |fl 1 1 |sty 1 e Iflxes... | | 

j-ldt; 

•f ul=u2fn - 1) 
ui 1=27 

' U2=|ul(n - 1!) 
ui2=21 
u3=| 
ui3= 
u4= 
ui4= 
u5= 

-u2(n- 1)1 

u3Cn>= 

MAIN RAD 

EKACT SEC: 





1 

h 

ul 

u2 




0. 

27. 

21. 



i. 

21. 

6. 




2. 

6. 

15. 




3. 

15. 

9. 




4. 

9. 

6. 




5. 

6. 

3. 




6. 

3. 

3. 




?. 

3. 

0. 




ln=0. 1 

IMAIN 

RAD EKACT 

SEC: 






H-fWl 

^p-4 1 

h 

u3 

u4 




0. 

1. 

.414214 



i. 

.414214 

. 171573 




2. 

. 171573 

.071068 




3. 

■071Q6S 

.029437 




4. 

. 029437 

.012193 




5. 

.012193 

.005051 




6. 

.005051 

. 002092 




?. 

. 002092 

8.67E-4 




1 

o 

II 

c 

IMAIN 

RAD EKACT 

SEC: 




friTmoT T F]: 
IZoonlEdit 

FH Y FEtY FbT Y F7 T 'I 

|h 1 1 Istyle Iflxes... | | 

jfLDTS 

ul=u2(.n - 1) 
ui 1=27 
u2=|ul(n - 1) 
ui2=21 
u3=u4(n - 1) 
ui3=l 

u4=u3(n - I) - 
ui4=J2 - 1 
u5=M 

-u2(n- 1)1 
2 u4(n - 1) 

u5Cn>= 

MAIN RAD 

EKACT SEC! 


Both examples from above performed with CellSheet Application of the Tl handheld: 


f FIT Y F£t Y F3t Y YFE Y Fb^ Y FFt T FB 1 

|F i le|Plot|Edit|:J: »Jo| $ |Funcs|stat|ReCalc| 

Q 

me a 

fl 

8 

C 

D 

E 

i 

27 

21 

1 

4X2)-1 


2 

21 

6 

4"<2>-l 

3-2*4" (2) 


3 

6 

15 

3-2*4" (2) 

5*4X2) -7 


4 

15 

9 

5*4X2) -7 

17-12*4X2) 


5 

9 

6 

17-12*4X2) 

29*4X2) -41 


6 

6 

3 

29*4" (2)-41 

99-70*4X2) 


7 

3 

3 

99-70*4X2) 

1 69*4" < 2) -239 


D7: =C6-2*D6 

MAIN RAD EKACT SEQ 


(Josef) 
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Tutorials for the NSpireCAS - Distributions 

Guido Herweyers, T3-Vlaanderen & T3-Wallonie, Belgium 
1. Sampling a normal distribution 


Press ((iT) to turn on the 
handheld. 

Press (fiT) for the home menu. 
Press <T> to open a new Lists and 
Spreadsheet page. 


Consider a population that has a 
normal distribution with mean 
/{ = 1 75 and standard deviation 
rr = 10 . 

We want to take a sample of this 
population with sample size 
n = 50 

Enter n in cell A1 and 50 in B1. 
Press for the context 

menu. 

Choose 7: Variables 
1 : Store Var 


Type n in the place of var. 

From now on, the value of this cell 
is linked to the variable n. 


Go to the gray position of column 
C (second line) to define the 
column and type: 

= randnorm(175,10,n) 

This will generate the sample in 
column C 


First you are asked to choose 
whether n is the name of a column 
or a variable. 

Press (J) and ▼ to choose 
variable reference, confirm with 
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The sample appears in column C 


Move to the top position to enter a 
name for column C, type 
sample 

Then press ▲ to select the whole 
column 


RAD BJTO REAL 


= 175.4895837!>!>S7 



B 


= 



RAD AUTO REAL 


E E G 


1 n 50 


Press Q|T)(^u) f or the context 
menu. 

Choose 9: Quick Graph 


RAD AUTO REAL 


.9: Quick Graph 


The window is split into two parts, 
by default a dot plot of the data 
appears in the data and statistics 
window. 


Press 

then 1 : Plot type 
3: Histogram 


Press 

then 2: Plot Properties 

2: Histogram properties 
1: Histogram scale 
3: Density 

We now have a histogram with 
density scale (total area = 1 ) 


c* oo 

Txl ffi 

cixrorooco' oooo 

155 165 175 185 
sample 


3: Histogram 


/ r 4: Normal Probability Plot 

[s*S: Scatter Plot 


5 q| | i Lire Plot 


o oo 

°ojS o° 

CEDO CJ>0CQi OQ3I! 

155 165 175 185 
sample 


ill, 1 : Histogram Seale > 
2: Bin Settings 


“•■-rr it'. 
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Press QjT) (S) for the tools menu, 
then 
5: Page 
1: Custom Split 


Press 4 for a wider histogram window 
and confirm with (j) 


Tools 


: Custom : 


RAD AUTO REAL 

OjO&oJ n . 


Press (™™), then 
3: Actions 
9: Plot Function 


The syntax for the normal density 
function with mean u = 1 75 and 
standard deviation o - 10 is 
normpdf(x, 175, 10) 


Confirm with 


Select the spreadsheet window with 

f cM 17 tab') 


Change the sample size into 100. 


| 157,. 

I 167 

I R*-- 

C I sample 


155 165 175 1S5 
sample 


1: Plot Type 
S Plot F ro pe rile : 
H 3: Actions j-gl: 
La, 4: Window ^2: 

X 


■ 



1 

n 

50 

2 



3 



4 



5 




, '.D AUTO REAL 
; >1 

Select All Points 
Add Movable Line 

Remove Selected 
Lock Intercept at Zero 
Regression 

Show Residual Squares 

show Normal PDF 
Plot Value 


KS9: Plot Function 


: Shade Under Function 

1 Start Over 


!!S 


n :=100 
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Go to the x-axis until +• appears, 
press (j) t0 get a closed hand £ 1 . 
Press ► or ^ to change the histogram 
width. 

Observe the influence of the histogram 
width. 


To change the window settings press 
then 

4: Window/Zoom 
1: Window settings 

choose Xmin= 145 and Xmax= 205 


Displace the function label. 

Select the spreadsheet window with 

( cW 17 tab ^ 


To take a new sample press Q^)(r) 

Take a few samples, the changing 
histograms illustrate the variability of the 
data. 


Change the sample size into 200, .... 
1000 . 

The histograms are getting closer and 
closer to the density function. 


Task: 

find out how you can produce the one 
variable statistics of the sample data in 
the spreadsheet window. 

The sample mean ~ and the sample 
standard deviation s are estimations for 
u - I 75 and rr - 10 . 

Observe: 

1 . The variability from sample to sample 
with the same sample size. 

2. The estimations are getting better 
with increasing sample size. 





1 1000 1 92.1 Title... One... 

156.. .. X ~1 175-.. 

180.. .. Ex 175... 
170.... Ex 3 3.07... 

185.. .. sx :... 9.75... 

E2 I =175.05767724365 


RAD AUT O REAL 

1 0.060 
^ 0.045 


® 0.030 


0.015 

0.000 


ormPd 



175 

sample 
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2 Binomial distributions 

Forty per cent of a sweet assortment are 'soft centres' and the remainder hard centres'. 

A handful of ten sweets may be regarded as a random sample from a huge pile of the 
sweets. 

a) What is the probability that there are three often soft centres' in a handful? 

b) What is the probability that there are at most four often 'soft centres’ in a handful? 


Press (^onp to turn on the handheld. 
Press (fiT) for the home menu. 

Press (J) to open a calculator page. 


Press then 
7: Statistics 
5: Distributions 
D: Biniomial Pdf 


The number of trials is n=10 
The success probability is p=0.4 
The X-value is 3 


If a random variable X has a binomial 
distribution with parameters n and p. then 

P(X=x) = binompdf(n.p.x) 

a) P(X=3) = 21.5% 

P(X<x) = binomcdf(n.p.x) 

(the cumulative distribution function) 



b) P(X<4) =63.3% 


2/90 
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3 Poisson distributions 


If a large grass lawn contains on average 1 weed per 600 cm 2 , what will be the distribution of 
the total number X of weeds in an area of 400 cm 2 ? 

a) Find P(X=2) 

b) Find P(X<3) 


A suitable model for this situation is a Poisson distribution whose mean is A = 


400 

600 


Press to turn on the handheld. 
Press (ffl) for the home menu. 

Press to open a calculator page. 


Press (™™), then 
7: Statistics 
5: Distributions 
H: Poisson Pdf 


ft Homy 


1 -ix 

L 

SI 

S 

E 

alculaton 

2:Graph$ & ... 

3:Liste & Sp... 


11 


m 

E 

4:Noie 

IS 

5;Dala ft Sta... 

6:Nri/» Doc... 


25 



<£> 

| 7 Docu... 

8:S/stem Inlo 

9:Hints 

1 Add a nei(V page with a Calculator application to 1 

| the open document. 

1 


ft/ 1- ArUr 


6: Inverse t 
7:x’ Pdf 
8: x’Cdf 

9: Inverse X s 

A:F Pdf 
B;F Cdf 
Converse F 
3: Biromlal Pdf 
E: Binomial Cdf 
F: Geometric Pdf 
G:Geometric Cdf 


a 


tat Calculations... 

► 

tat Results 


ist Math 

> 

ist Operations 

► 

■istrlbuCons... 

► 

ionlidence Internals... 

► 

tatTests... 

► 


A = 2/3 and X = 2 


PAD AUTO REAL 


A! | 2/3 

X value: 2 


OK 


Cancel 




If a random variable X has a Poisson 
distribution with mean A , then 

P(X=x) = poisspdf(i ,x) 

a) P(X=2) = 1 1 .4 % 

P(X<x) = poisscdf(/l ,x) 

(the cumulative distribution function) 

b) P(X<3) = 99.5% 


u i..- 


RAD AUTO REAL 


poissFdf ~ .2 
3 I 


.114093 


poissCdf | — ,3 


.995142 


You can find much more materials at http://www.t3vlaanderen.be and http://www.t3wallonie.be . 
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What's the Time, Grandie? 

Josef Bohm, Wiirmla, Austria 

Maybe that some of you will remember DNL#45 where I introduced proudly our first granddaughter 
Kim. Later I presented a picture with Kim and the TI-89 (with Grandpa's eyeglasses, of course). Now 
Kim is attending the grammar school in Tulin, Lower Austria, and sometimes - it is not too often - she 
asks her grandpa for support in mathematics. 

Below is Kim together with her sister Yvonne 
sitting on my desk some years ago. 

At the right you can see Kim and Yvonne (both 
left) together with their younger sisters and her 
mother Astrid (our daughter). 




It was last year when Kim needed some exercising for calculations with times. I took the occasion to 
extent my collection of training programs and wrote a DERIVE utility for adding and subtracting times 
and for adding times to certain day times in the week. 

•y = ■+ « % x |start= 

D £ H 5 Ej Ej M [::!] = » s Su B lim a j 2 n “p k & I 

#1: LOAD(D: \DOKUS\SCHULE\ski 1 ls\ti me US . mth) 

Trainer for calculating times 

ti(n)= gives n problems to add and subtract times and 
time_ex_res= shows the correct answers. 

time= gives a day of the week together with the time of the day 
and a time to add or to subtract. 
res= gives the correct answer. 

time_ex(n)= offers n examples in form of a table and 
time_ex_res= shows the correct answers. 
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All my training programs show the same structure (except one which treats exercising with Venn- 
diagrams for visualising set operations): 

The files are MTH- files and should be loaded as Utility- files. Then simplify the command start. You will 
be presented the instructions how to use the file (see page 15). The “trick” is to use the DISPLAY com- 
mand for this purpose (see the code of start_() on the next page). start_() contains the simplification of 
dummy:=random(0) which makes sure that we will be offered new problems at every run of the utility. 

A session could - according to the instructions given above - start and run as follows: 

#2: start = 

5 hrs 53 min IS sec + 13 hrs 37 min IS sec 
#3: ti(3) = 17 hrs 26 min 11 sec - 4 hrs 35 min 20 sec 

20 hrs 9 min 39 sec + S hrs 17 min 32 sec 

19 hrs 30 min 36 sec 
#4: time_ex_res = 12 hrs 50 min 51 sec 

. 2S hrs 27 min 11 sec 

#5: time = Fri 4 hrs 30 min am - 7 hrs 7 min 

#6: res = Thu 9 hrs 23 min pm 

#7: time = Thu 4 hrs 46 min pm + 17 hrs 38 min 

#8: res = Fri 10 hrs 24 min am 

Thu 12 hrs 23 min am + 16 hrs 0 min 
Thu S hrs 35 min pm - 5 hrs 44 min 

#9: time_ex(5) = Mon S hrs 53 min am - S hrs 55 min 

Wed 3 hrs 49 min am + 15 hrs 12 min 
_ Mon 3 hrs 55 min am + 6 hrs 24 min 
Fri 4 hrs 23 min am 
Thu 2 hrs 51 min pm 
#10: time_ex_res = Sun 11 hrs 5S min pm 

Wed 7 hrs 1 min pm 
. Mon 10 hrs 19 min am 

Kim learns English - and she likes English more than mathematics - so I could not resist to include 
the “am - pm” notation of the times of the day. Additionally it was nice to program this feature. 


The next page shows the utility file: 
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#1: 

[days 

:= [Mon 

i Tue , 

Wed 

, Thu, 

Fri, Sat, Sun]] 

#2: 

[time_ 

.res :=, 

t_res 


time_ex 

_res :=] 

#3: 

I 1 

r+ 

H 1 

i i 

hi :=, 

ml :=, 

h2 

:=, m2 := 

1 1 

1 [ 

CL 

□ 


start_(dummy) 

Prog 

dummy := RANDOM (0) 

DISPLAY("Trainer for calculating times") 

DISPLAYC" 1 ) 

DISPLAY("ti(n)= gives n problems to add and subtract times and") 
DISPLAY("time_ex_res= shows the correct answers.") 

#4: DISPLAYC" 1 ) 

DISPLAYC 1 time= gives a day of the week together with the time of the day") 
DISPLAY("and a time to add or to subtract.") 

DISPLAY(" res= gives the correct answer.") 

DISPLAYC" 1 ) 

DISPLAY("time_ex (n)= offers n examples in form of a table and") 
DISPLAY("time_ex_res= shows the correct answers.") 

DISPLAYC 1 ") 

#5 : start start_() 

zeit_(days, timel, time2, hl_) := 

Prog 

tl := RANDOM (7) + 1 
hi := RANDOM (24) 

hl_ := IF (hi >12, hi - 12, hi) 
sp := IF(hl > 12, "pm", "am") 

#6: ml := RANDOM (60) 

timel := APPEND(daysj.tl, " ", STRING(hl_), " hrs ", STRING (ml), " min ", sp) 
h2 := RANDOM (24) 
m2 := RANDOM (60) 

time2 := APPEND(STRING(h2) , " hrs ", STRING(m2) , " min") 
op := [" + ", " - "]j. (RANDOM (2) + 1) 

APPEND(timel , op, time2) 

#7: time := zeit_() 


res_(st 

d, erg, mr, hr, 

tr 

^) 

Prog 




st : 

= tl-24-60 + hi- 60 

+ 

ml 

d := 

h2 -60 + m2 



erg 

:= IF (op = " + ", 

st 

+ d, 

rnr : 

= M0D(erg, 60) 



erg 

:= (erg - mr)/60 



hr : 

= M0D(erg , 24) 



erg 

:= M0D((erg - hr)/24 

7) 

sp : 

= IF(hr > 12, "pm" 

, 

’am") 

erg 

:= IF(erg = 0, 7, 

er 

3) 

hr : 

= IF(hr > 12, hr - 

12, hr 


APPEND(days|erg , " ", STRING(hr), " hrs ", STRING(mr), " min ", sp) 
#9: res := res_() 

time_ex(n, n_, table) := 

Prog 
n_ := 1 
table := [] 
time_ex_res := [] 

#10: Loop 

If n_ > n exit 

table := APPEND(table , [[time]]) 
time_ex_res := APPEND(time_ex_res , [[res]]) 
n_ : + 1 
RETURN table 
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The next expression is expression #11.1 removed the expression number to have a better screen shot 
of the extended function. 

time_(secs, min, hrs, timl, tim2, tim3, signs, si, tsk, pn) 

Prog 
Loop 
pr := "" 

signs := VECTOR (1 - 2 -RANDOM (2), k, RANDOM (2) + 1) 

tsk := VECTOR (RANDOM (60 -60-24) + 1, k, DIM (signs) + 1) 

signs := APPEND([1], signs) 

secs := signs -tsk 

t_res := secs 

timl M0D(secs, 60) 

t_res := (t_res - timl)/60 

tim2 := M0D(t_res, 60) 

tim3 (t_res - tim2)/60 

t_res := APPEND(STRING(tim3) , n hrs n , STRING (tim2) , 11 min ”, STRING (timl) , n sec 1 ’) 

If secs > 0 
Loop 

If DIM (tsk) = 0 

RETURN REST (REST (REST (pr))) 
secs ;= M0D(tsk|l , 60) 
tskj.1 := (tskrl - secs)/60 
min := M0D(tskj,l, 60) 
hrs := (tskj.1 - min)/60 
si := IF(signs|l = 1, " + 11 , " - n ) 

pr := APPEND(pr, si, STRING(hrs) , 11 hrs 11 , STRING (min) , 11 min ", STRING(secs) , " sec ") 
tsk := REST (tsk) 
signs := REST(signs) 

ti(n, n_, table) := 

Prog 
n_ := 1 
table := [] 
time_ex_res := [] 

#12: Loop 

If n_ > n exit 

table := APPEND(table , [ [time_() ] ]) 
time_ex_res := APPEND(time_ex_res , [[t_res]]) 
n_ : + 1 
RETURN table 


Among the files accompanying this DNL you can find zeit . mth, which provides a German version 
without the “am-pm”-notation. It should be no problem for non-German users to adapt the file for their 
own language if necessary, (tage : = [Mo, Di,Mi, Do, Fr, Sa, So] must be changed!) 


#12: zeitb(5) = 


#13: zeitb_res = 


Do 8h Onrin + 3h Omin 
So 2h 2min - 12h 59min 
Do 9h 5min + 15h 18min 
Fr 20h 43min + 12h 13min 
Do lOh 29min + llh 2min . 
Do llh Omin 
Sa 13h 3min 
Fr Oh 23min 
Sa 8h 56min 


. Do 21h 31min . 
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An Interesting Property of a Triangle 

Peter Liike-Rosendahl, Germany 

In 1968 the following theorem was presented as a problem in the journal American Mathematical 
Monthly. 

Given is a triangle ABC. We raise squares CBED, ACFG and BAHK over the sides 
of the triangle outwards. Then we draw the parallelograms FCDQ and EBKP. Tri- 
angle PAQ seems to be right and isosceles as well. Prove this! 

(It might be nice to pose the problem to the students changing the last sentence: What can 
you say about triangle PAQ ? Prove your conjecture? Josef) 


The first step (finding a conjecture) is a typical 
task for working with a dynamic geometry pro- 
gram. 

We can do it on the PC using TI-NspireCAS or 
Cabri or GeoGebra or any other DGS. 

But we can do it on the handheld, too using either 
the TI-92 / Voyage 200 with their Cabri- 
Application or the TI-Nspire handheld device. 

The next step could be verifying the conjecture 
using special data for three points and verifying 
the fact that the other two possibilities drawing 
the triangle will give the same result. 



I will proceed with DERIVE using the slider bars keeping the coordinates of the triangle vertices gen- 
eral and then continue proving the two properties. 

#1: CaseMode := Sensitive 

#2: InputMode := Word 

#3: triangle(U, V, W) := [U, V, W, U] 

#4: [A := [0, 0], B := [xb, 0], C := [xc, yc]] 

Introduce sliders forxb, xc and yc from -10 ... +10 

#5: triang1e(A, B, C) 


Plot the triangle 
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-3 -2 -1 

perp(U, V) 


10 


# 6 : 


(V - U) , - CV - U) 

2 1 . 


square(U.V) = square over UV outside of the triangle (neg. direction of rotation) 
#7: square (U, V) := [U, V, V + perp(V, U) , U + perp(V, U) , U] 

#8: [square(A, C) , square(C, B) , square(B, A)] 

Plot the squares 
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Ppt(U,V,W) = point of parallelogram with point of triangle W as vertex 
#9: Ppt(U, V, W) := W + perp(W, U) + perp(V, W) 

par(U,V,W) = parallelogram with vertex W 

#10: par(U, V, W) := [W, W + perp(V, W) , Ppt(U, V, W) , W + perp(W, U) , W] 

Plot the parallelograms 



#11: [par (A, B, C) , par(C, A, B), par(B, C, A)] 

#12: goaltr (U , V, W) := [U, Ppt(U, V, W) , Ppt(W, U, V), U] 
#13: goaltr (A, B, C) 

Plot one of the final triangles 

#14: P0LYG0N_FILL([A, B, C, A]) 

Plot to fill the initial triangle 

#15: [goaltr(B, C, A), goaltr(C, A, B)] 
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#15: [goaltr(B, C, A), goaltr(C, A, B)] 

Plot the other two final triangles 



Prove the properties of the final triangles: 

#16: (PptCA, B, C) - A) ■ (Ppt(C , A, B) - A) = 0 

#17: | PptCA, B, C) - A| - | Ppt(C , A, B) - A| = 0 

#1S: (Ppt(B , C, A) - B). (Ppt(A, B, C) - B) = 0 

#19: | Ppt(B , C, A) - B| - | PptCA, B, C) - B| = 0 

#20: |Ppt(C, A, B) - C)-(Ppt(B, C, A) - C) = 0 

Ppt(C , A, B) - C| - |Ppt(B, C, A) - C 


# 21 : 


0 
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I asked myself if this special triangle would have a nice result for its area and reproduced the 
construction using GeometryExpressions. Josef 



At the bottom you can see a part of the formula for the area with given sides a, b and c of the 
initial triangle. It is a “nice” formula, indeed. I can export this formula to DERIVE, 
MATHEMATICA, ... - and with the new GE-version to TI-NspireCAS, too. 

JU>s[ 


((Abs[ ((a) A (2) ♦ ((b) A (2). (-!>) + (c) A (2))]) A ((-!)) * 


(Sin[ 


< (ArcTan[((((a) A (2) * ((b) A (2) * (-1)) + (c) A (2))) A ((-1)) * (((a« (-1)) + b + c)) A (1/2) * ((a + (b + (-1)) + c)) A (l/2) + 

((a + b + (c.(-l)))) A (l/2).((a + b + c)) A (l/2),(-l))]„(-l» + 

RrcTan[ ( ( ( ( (a) A (2) * (-1) ) + (b) A (2) ♦ (c) A (2))) A (<-l>) , (((a, (-1)) * b * c)) A (1/ 2) » ((a + (b , (-1)) + c)> A (1/ 2) , 

((a+ b + (c * (-l)))) A (l/2) * ((a + b + c)) “(1/2))])]) A ((-l)) » 

<((a) A (6). (-1)) ♦ ((b) A (2) * (a) A (4) * 3) + ((b) A (4) * (a) A (2) * (-3)) + (b) A (6) + ((c) A (2) * (a) A (4)) + 
l((c) A (2), (b) A (2), (a) A <2), (-2)) + «c) A (2) » (b) A (4)) + (<c) A (4) „ (a) A (2)) + ((c) A <4). (b) A (2) , (-1)) * 

((c) A (6), (-1)) * ((((a,(-l)) + b+ c)) A (l/2)*((a + (b * (-1)) + c)) A (1/ 2) , ((a + b . (c , (-1)))) A (1/ 2) , 

((a + b + c)) A (1/ 2) * (a) A (4) it (- 1)) + 

((((a* (-1)) + b + c)) A (l/2) * ((a + (b , (- 1)) + c)) A <1/ 2) , ((a * b * (c , (-1)))) A (1/ 2) , ((a + b * c>) A (1/ 2) , (b) A (2) , (a) A (2) , 2) * 

((((a.(-l)).b.c» A (l/2),((a* (b.(-l)).c))‘(l/2).((a.b. (c . (-1)))) A (1/ 2) „ ((a +b + c)) A (1/ 2) * (b) A (4) „ (-1)) + 

((((** (-1)) + b + c)) A (l/2) * ((a + (b„(-l))+c)) A (1/2). ((a + b. (c . (-1)))) A (1/ 2) . ((a +b . c)) A (1/ 2) » (c) "(4)) ♦ 

(M>s[((a) A (2) + ((b) A (2).(-l)) + (c) A (2))]»(((a,(-l)) tb ,c»«(l/2) .((a. (b . (-1)) ♦ c» A (1/ 2) , 

((a,b*(c,(-l)))) A (l/2),((a*b.c)) A (l/2),(a) A (2)) + 

(»>s[((a) A (2) + ((b) A (2).(-l)) + (c) A (2))]»(((a, (-1)) * b + c)) A (1/ 2) , ((a* (b . (-1)) ♦ c)) A (1/ 2) . 

((a+b + (c . (~l)))) A (l/2) . ((a + b + c)) A (l/2) * (b) A (2) . (-1)) * 

(Sbs[((a) A (2) * ((b) A (2).(-l)) 1 (c)'(2))].(((a*(-l))*b*c))*(l/2).((a, (b . (-1)) ♦ c)) A (1/ 2) . 

((a + b + (c * (-l)))) A (l/2) * ((a + b + c)) A (l/2) . (c) A (2) * (-1)) + 

(Bbs[(((a) A (2). (-1)) + (b) A (2) + (c) A (2))] * 

(((a) A (4) * (-1)) + ((b) A (2).(a) A (2).2) + ((b) A (4) . (-1)) + (c) A (4) + 

((((a.(-l)) + b+c)) A (l/2).((a + (b.(-l))+c)) A (l/2).((a + b + (c.(-l)))) A (l/2).((a+b+c)) A (l/2).(a) A (2)) + 
((((a.(-l)) + b+c)) A (l/2).((a + (b . (-1)) + c)) A (1/ 2) . ((a + b + (c . (-1)))) A (1/ 2) . ((a +b + c)) A (1/ 2) . (b) A (2) . (-1)) + 
((((a.(-l)) + b+ c)) A (l/2)*((a + (b.(-l))+c)) A (l/2).((a + b + (c.(-l)))) A (l/2).((a + b+c)) A (l/2).(c) A (2)) + 
(Hbs[((a) A (2) + ((b) A (2).(-l)) + (c) A (2))].(a) A (2).2) ♦ 

(Bbs [((a) A (2) ♦ ((b) A (2).(-l)) + (c) A (2))].(((a.(-l)) + b+c)) A (l/2) . ((a + (b . (-1)) ♦ c)) A (1/ 2) . 

((a+b + (c.(-l)))) A (l/2)*((a + b + c)) A (l/2))))). 

((a) A (4) . ((b) A (2) . (a) A (2) . (-2)) + (b) A (4) ♦ ((c) A (2) . (a) A (2) . (-2)) ♦ ((c) A (2) . (b> A (2) , (-2)) ♦ (c) A (4)) . 

(((a) A (2) + ((b) A (2). (-1))+ ((c) A (2).(-l)))) A ((-l)). (((a.(-l))+b+c)) A ((-l/2)) . ((a+(b. (-1)) + c)) A ((-1/ 2)) , 

((a+b + (c * (-!))))"( (-1/2)) * ((a+ b + c)) A ((-1/ 2)) * (b) A ((-!)) * (a) A ( (— 1) ) +1/8)] 


The MATHEMATICA - output of this formula does not improve its appearance! 
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Peter delivered a formula for the area: 


[A := [0, 0] , B := [b, 0] , C := [c, d]] 

Ppt(A, B, C) = [c, b + d] 

Ppt(B , C, A) = [-d, c - b] 

Ppt(C, A, B) = [b + d, -c] 

Area of the initial Triangle 
1 

CR0SS([b , 0] - [0, 0], [c, d] - [0, 0]) 

2 

Area of the goal triangle 


I b ■ d | 
2 


-CROSSCtc, b + d] - [0, 0], [b + d, -c] - [0, 0]) 

2 


2 2 2 
b + 2‘b-d + c + d 


2 


or 


1 

-CROSSCfc, b + d] - [b + d, -c], [b + d, -c] - [0, 0]) 


2 


2 2 2 
b + 2*b*d + c + d 


2 


I repeated the construction using GeometryExpressions with generalized coordinates of the 
edges A, B and C and hoped for a nicer result: 




Derive Input 
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Then I exported this result to DERIVE (other possibilities are to MAPLE, MATHEMATICA, 
Tl-Nspire, Maxima, Mupad). 


#1: InputMode := Wbrd 


# 2 : 


2 2 

2 xl -1 x2 -1 2 

xO + xl.xO.(-l) + + x2-x0.(-l) + + yO-xl-(-l) + y0-x2 + yO + yl-xO + yl.x2.(-l) + yl.yO.(-l) 

2 2 


2 2 
yl -1 y2 -1 

+ y2*x0'(-l) + y2-xl + y2-y0*(-l) + 


2 2 2 2 
2 xl x2 2 yl y2 

#3: xO - xO-(xl + x2 - yl + y2) + + xl*(y2 - yO) + + x2*(y0 - yl) + yO - yO*(yl + y2) + + 

2 2 2 2 


2 2 2 
2 b -1 c -1 2 0-1 

0 + b-O.(-l) + + c ■ 0 ■ (-1) + + 0 ■ b • (-1) + 0-c + 0 + 0-0 + O-c-C-l) + 0-0. (-1) + + d-O.(-l) + 

2 2 2 


d ■ b + d ■ 0 ■ (-1) + 


2 

d -1 


2 



I substituted Peter's special values and could confirm Peter's results. This formula looks 
much better, isn't it? 

Many thanks to Peter for his inspiring contribution. 

By the way, I gave this problem to my students at the Technical University Vienna and sur- 
prisingly many of them tried the proof not using vector calculations but worked with trig func- 
tions (sine and cosine rule ...). 

Website of GeometryExpressions: 

www . geometryexpressions . com 


Everyone knows that mathematics offers an excellent opportunity to 
learn demonstrative reasoning, but I contend also that there is no 
other subject in the usual curricula of the schools that affords a com- 
parable opportunity to learn plausible reasoning. I address my self to 
all interested students of mathematics of all grades and I say: Cer- 
tainly, let us learn proving, but also let us learn guessing. 


G. Polya (1954). Mathematics and Plausible 
Reasoning. Princeton, NJ. Princeton Univer- 
sity Press. 
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Differential Equations Made Easy 

Review and Description of some features, Josef Bohm 


DUG Member Nils Hahnfeld gathered a bundle of programs and functions into one application 
DEQME v.7.0, which can be downloaded and purchased from www.ti89.com . The application is 
available for TI-89, TI-92+ and Voyage 200. 

Nils and I had an extended exchange of emails and files. In some details he could refer to earlier DNL 
contributions. 


I will give an overview about the many features of DEQME. 


First of all see screen shots which show all options offered in the various menus reaching from Basics 
of 1 st order DEs over PDEs to special DEs, Laplace Transforms and Eigenvalues. 




FEt T F3t 
2. Order 1 0 r d e r n 

FHtT FSt T FbT "\t? 'I 

PDE|Transf |More|> | 

fltN<x,yJy'+MCx,y)=0 
B:l. Order DE Checker 
C:RL Circuits 
D: Bernoul 1 i 
E:Clairaut Equation 
F: Lagrange Equation 
G: Linear Fractions 
H: Slope Field 
I: Par tic. Solution 
J: Euler's Method 
K: Runge-Kutta 

is 





TYPE OF; USE 4*t+ + [ENTEM=UK AND [ESC]=CANCEL 


r fit 

|1 . Order 

F3t T FHtY FEt Y FbT "if? ] 

^!5®HS0rder n |PDE |T ransf |More | ► | 


lsflny 2. Order DE 


Differ 

2:Honogeneous 
3 ! N o n - H o n o gene o u s 
4: 2. Order DE Checker 
5: IUP So Iyer! xl*x2 
b : U a r i a i i o n o + 1 Pa r a . 
7 : Un determ i ned Coeff 
S ! Be sse 1 Eq u a t i on 
9 ! Re d u c t i o n of 0 r d e r 
fi: Legendre DEQ 
E: : C a u c h y - Eu 1 e r D E Q 



TYPE DF; USE 4»ti * [ENTER]=DK AND [E;C]=CHHCEL 



Differential 


2:2x2 Linear 
3: X 1 =R*X 
4: X 1 =fl*X + F 
5: Separable PE 
Made Easy 


TYPE DF; USE 4»ti * EENTER]=DK AND [E;C]=CHHCEL 


( iTF Y Y F3 t ' 

|l ■ Order |2. Order Idrder n 

Differential Eqi 
Made Easy . 



2: LaPlace Equation 
3: Diffusion Equation 
4: Wave Equation 
5:Helnholtz Equation 
6 : Po i sson Equation 


TYPE OP USE 4»ti * EENTER]=DK AND [E;C]=EHMCEL 


f Fi t Y f£t T F3t \^m 
\ 1 . Order |2 . Order |0rder n PDE iff 

FfiT "it? -i 

efflaaMnre ► 




Differential Eqr 
Made Easy 

2: Solve DE 
3: Laplace 
4: Inverse 
5: Inv. Lac 
6: FourierE 
7: Partial 

.0 Laplace 
Transf orris 
Laplace 
>lace-STEP 
jeries-STEP 
Fractions 




TYPE □ F; USE + [ENTEF;]=DK AND [ESC]=CANCEL 


r Fi T j F Y F3 t T FHt 

1 1 . Order |2 . Order |0rder n |PDE 

FSt 'll 

T ransf u 


hr? i 
1 






Differential Equ; 
Made Easy 

2: Rate Problens 
3: Salt in Tank 
4: Wronskian 
5: Flnnihi lator 
6: Picard's Method 
7: Error Function 
8: Phase Line 






TYPE OF; USE + [ENTEF;]=UK AND [ESC]=CANCEL 
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I open the first menu point under FI 1 .Order, 1 : Basics and want to inform about General Sol(ution): 



Then I switch back to FI , 2:Any 1 .Order DE and would like to solve the differential equation 
2y - x - 5 


y = 


2x - y + 4 


; y(0) = 4 . 



I don't receive the solution because the TI built in desolve cannot solve this kind of DE. Under FI 1 
can find the option G:Linear Fractions. So I try again: 



This is the general solution, but how to obtain the special solution? The solution which is presented in 
the Prgm 10 screen is stored as res and can be recalled in the Home screen. In F7 Menu you can find 
the respective note. 
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11 


J| 


Just for checking the result I load DERIVE and apply the L I N_FRAC -function: 


LIN_FRAC 


2-y - x - 5 
2-x - y + 4 


-1, 2, -5, 2, -1, 4, x, y, 0, 4 


LN 


x - y + 3 


x + 1 


3- LN 


x + y - 1 


x + 1 


= LN(x + 1) 


3 - LN (35 


2 2 2 

After some manipulations I obtain the same solution! Well done, DEQME!! 

3 2 2 2 3 2 

#11: x + 3-x -y - 3-x + 3-x-y - 6-x-y + 30- x + y - 3-y - 24-y 


-SO 


In case of non intersecting linear functions in numerator and denominator DERIVE provides a special 
utility function FUN_LIN_CFF_GEN. Nils implemented this special case as you can see below: 



Pi 

□ 





f Linear Fractions 


DE: : |h'=<2x+li-2)^C6x+3u-10) 1 

Exl: y '=(2x+y-3:>/(x-y+2:> 

CEnter=0K CESC=CRNCEO 




MAIN RAD AUTO FUNC 1/30 


I try solving another DE using the first menu point “Any 1 . Order DE” 1 + y 2 + x ■ y ■ y' = 0; y( 3) = 2. 



Enter any 1. Order DE. 

□ 



DE: : ll+y^+x+y+y 1 =0 1 


Exl : s ' =s^2 
Ex2: x*y'=y 

I ride p. var.=! |x I 

Dep. uar.=: |y i 

OPTIONAL: Initial condition 

at x or t or . . = : 13 i 

y or s or . . = • |2 1 

t CEnter=0K <fESC=CANCELV 

ImAIN RAD hUTD FUNC 0/30 |i 



And I’d like to plot its slope field and the special solution. 


Pi 

□ 

1 

Enter DE to Graph slope field. 


dyxdx=: 1 -a+y''2)/Cx*y>l 1 

N0TE1: Use variables x and y. 
N0TE2: Only dyxdx on left side. 
Exl: dyxdx=x*y 
Ex2: dyxdx=x-y 

CEnter=0K <TESC=CANCEL^ 




MAIN RAD HUTD DE 1/30 


’ FIT 

1 . Order kwJ 1=1 J L«J Li Ws 



IS 

■esd uitiai iiiuuia. r 




MAIN RAD AUTO FUNC 1/30 


As you can read in the bottom line, the GRAPH Mode changed from FUNC to DE. 
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!■"- '"iirY 

? ] • Y : : v Y ■ 

y: 

# :■ Y #: i'T ’■ 

□ 

Li 

f Port. Solution - Entir PE ond point \ 


1 dyxdx=: l-Cl +u''2V 

<x*y) || 



N0TE1: Use uariables 
Ex: dyxdx=x*y 

x and y. 



X= 5 

13 


1 



Ex: 

x=0 





y=: 

la 


1 



Ex: 

y=2 





CEnter=OK 0 


("ESC=CflHCEL">j 







IMhIN 

RAD AUTO 


[iE i/30 


F£T Y F3t 

2 . 0 r d e r 1 0 r d e r n 

FHtY FFt T FhT 'it? i 

PDE|Transf |Morek | 

FltN(x, ypy 1 +M(x, y )=0 
B:l. Order DE Checker 
C:RL Circuits 
D: Eernoul 1 i 
E:Clairaut Equation 
F: Lagrange Equation 
6: mnear^_^ract ions 

is 

V3ME?HI 1 V L 


R: Runge-Kutta 
Ls Midpoint- 



MAIN RAD AUTO DE 4/30 


fl^lz^nlTrlcelRe^aphlN^hlD^llVIlC 1 

ss-/- / / 1 \ 

; y 7 

— '■r' /': a / ./ 

T\V \ \ *W n. ^ 
\ \ \ \ 

V 

^ ^ \ -\ 
\\\ 

- -7 \ \ \ .l 

-- > - .V.\. \ V 

//a > 

/. / / a a.-a-a-a-tt- 

\. i. / / y.v. ^ ^ 

1 MAIN RAD AUTO DE 1 


3x v 

y' = — , v(-2) = 4. Solve and plot the solution together with the direction field. 

2x +y 



Enl4F OtiV i. OrdiF DE. 

□ 



DE: : l3x*y-C2x''2+y )*y ' =0 1 


Exl : s'=s A 2 
Ex2: x*y'=y 

Indep. uar.=: lx 1 

Dep. uar.=s |y 1 

OPTIONAL: Initial condition 

at x or t or . . = : 1 "2 1 

y or s or . . = : |4 1 

k CEnter=OK 5 CESC^FlNCEDj 

ImAIN RAD AUTO FUNC fi/30 t 



fl^Iz^lTr^celReg^aphlM^hlD^JlVllC 1 

vx'. \ VV \ A 

^ ^ ^ \ S ^ ^ \ 

■— — - ■— ^ \ \ \V 

v \ ^ 

/■ //? / A ^ ^ X 1 

ff'/y/ssss's' 
y/ / ,x> 

7 

< A A. A / / 7 A 

s.s / / S./.f. i .1 A 

/■ \ 

y \\\ 

I.\\ \ >>:•-• 

/. 1. t \ \ V .v 

IMAIN RAD AUTO DE 


H i (A T : ii- T =-^ r FE 1 ' 

□ 

Li 

f Fart. Solution - Enter DE ond point Cx^p). ^ 


duxdx=: l3x*yx(2x''-2+y:J 1 


NOTEls Use variables x and y. 
Ex! dyxdx=x*y 

x=: 1-2 1 

Ex: x=0 

y=: [4 1 

Ex: y=2 

, ("Enter =0K "> CESC=CflNCEL"> J 




IMAIN RAD AUTO DE fi/30 1! 



Runge-Kutta is implemented in order to find 
numerical solutions. 




v Y i i V Y : Y * i‘ Y # :i'r 

P 0 


Run34-Kutto M it. hod: 


Flpproxinates y coordinates 
using: 

k_l=h*f (jf(n),u(n)) 
k_2=h*f(x(n>+hx2,y(n>+k_l/2> 
k_3=h*f(x(n)+hx2,y(n)+k_2/2> 
k_4=h*f (x(n!>+h, y (n^ + k_3J 
Then y(n+l>=y(n> 

+ ( Ll+2*k.2+2+k.3+k.4)/6 

("Enter=0K "> ("ESC=CflNCEL"> ^ 





MAIN 

RAD AUTO FUNC 2/3 0 



«^H|" F2t Y F3t 

mjhi[=jm 2 . 0 r d e r 1 0 r d 0 r n 

f FHtY FFt Y FhT "if? A 

PDE|Transf | Nor -0 | 

EtN(x, y!>y 1 +M(x, y I>=0 
E:l. Order DE Checker 
C:RL Circuits 
D: Eernoul 1 i 
E:Clairaut Equation 
F: Lagrange Equation 
0: Linear Fractions 
H: Slope Field 
I:Partic. Solution 
J: Euler^s Methoc^^^ 

is 

L ! fUdpoTnt 




MAIN RAD AUTO FUNC 2/30 
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I'd like to see the calculation steps: 
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X V 

y' = — ; v(l) = 2. This is a homogeneous DE. 

y-x 

Again I'd like to I follow the steps. Then I will compare with deSolve. 




How to find an integrating factor— if there is any? DEQME gives the answer: 



Usinq Initial 

conditions: 

G = 2- 
Thus, 

a + 2^3 

P a r t i c u 1 a r so 1 u i i o n : 


_.y _ 

y3 -2a + > 

2 

Done. 

a x X 

—r ~ jl. -3 T . r 

3-x J 
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Titbits(37) - Factoring Integers with DERIVE 2 

(c) Johann Wiesenbauer, Vienna University of Technology 


This is a sequel of my last titbits where I started with the discussion of several algo- 
rithms for the factorization of integers, namely Pollard's two factoring methods (rho 
and p-1) as well as Lenstra's ECM. In the following I'll talk about another basic idea 
when it comes to factoring integers, which deals with nontrivial solutions x,y of the con- 
gruence x~2 = y~2 mod n, where n is again the positive integer to be factored and "non- 
trivial" means here that x t ±y mod n. Since due to these conditions n is a divisor of the 
product (x+y)(x-y), but not of its factors x+y and x-y, it easily follows from this that 
both gcd(x+y,n) and gcd(x-y,n) will be nontrivial divisors of n. 


This simple idea can be exploited in several ways, but in the following I'll focus on the 
so-called quadratic sieve only. Here one considers polynomials Q(x) in Z[x] of the form 


Q(x) = (ax+b) A 2-n = af(x), with f(x)=ax~2+2b x+c and c=(b~2-n)/a an integer 

and their values in the interval I=[-m,m] that contains both real zeros and has the prop- 
erty that the extrema of Q(x), located at the endpoints of I and near its center, have 
about the same absolute size. This will guarantee that the values of f(x) for x in I are 
relatively small, which is important as we want f(x) to be B-smooth for a bound B > 0 of 
moderate size, that is we are only interested in values of x in I such that p <=B for all 
prime divisors of f(x). Well, actually we also consider also the cases where this is "al- 
most" true, in the sense thta it is true except for at most one not too divisor p, exceed- 
ing B not too much, say B< p <=10*B. The underlying idea is that if we find two represen- 
tations zVZ - yl * p mod n and z2 = y2 * p mod n with yl and y2 both B-smooth numbers 
and some common prime p > B, then combining them like this 


(zl*z2/p) / '2 = yl*y2 mod p 


yields a representation where now the right side is perfectly B-smooth. This situation 
will occur surprisingly often, which is also known as "birthday paradox". 


As for the choice of a, there are several possibilities, but we choose here one of the 
simplest, where a = q~2 ~ sqrt(2n)/m for some prime q. Furthermore, we choose b to be 
a solution of x~2=n mod a with |b| < a/2. These choices of a will guarantee that the 
coefficient c of f(x) is an integer indeed, and that the extremal values of f(x), located 
at the endpoints and the center of I, have about the same absolute size. 
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Since Q(x) is a square mod n and a fortiori also mod p for any prime divisor p of n, we 
need only consider primes p for which n is a square mod p. This is always fulfilled for 
p=2 and for an odd prime p equivalent to jacobi(n,p)=l, where jacobi() is the so-called 
Jacobi symbol, which exists as a library function in DERIVE. Hence, if we consider the 
set 


FB = {-1,2} u {p<=B | p is an odd prime and jacobi(p,n)=l} 


then we want to factor f(x) using only numbers in FB. For this reason, FB is also called a 
factorbase w.r.t B, denoted by f in the program. 


Lets start our implementation with the fundamental routine setup(v,P) that sets the 
values of a number of important global variables, namely 


n = the number v to be factored, 

b = upper bound for the primes in f (if this parameter is omitted, it will be chosen by 
the program in a reasonable way) 

f = factorbasis consisting of -1,2 and the odd primes p below b such that n is a quadratic 
residue mod p, 

m = ceiling(b/2), which is chosen in such a way that the interval length of I=[-m,m] is 
about b, 

fp = product of all primes in f, which will be needed later for technical reasons 


User: 

setup(v, 3 := 0) := 

Prog 
n := v 

b := IF(p > 0, 3, CEILING(EXP(0.43-JCLN(n)-LN(LN(n)))))) 

# 1 : f := APPEND( [-1 , 2], SELECT(PRIME(q_) a JAC0BI(n, q_) = 1, q_, 3, b)) 

m := CEILING(0 . 3 • b) 
fp := n(REST(f)) 
true 

User=Si mp(User) : 

#2: setup(24961) = true 

User=Si mp(User) : 

#3: [n, b, f, m, fp] = [24961, 9, [-1, 2, 3, 5], 3, 30] 

The next routine getrel() is supposed to collect a sufficiently large number k of rela- 
tions, i.e. pairs (u i,v i), such that 


(u_iT2 = vj mod n, i=l,2 k 
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and vj is b-smooth, where u_i and v_i are obtained by using polynomials (ax+b^-n for 
various values of a and b as outlined above. The output of getrel() is a vector x_rel, con- 
taining the u_i , and a matrix E_rel, whose row vectors represent the exponents of the 
prime factorizations of the vj over the factor basis f. 


getre"l(k := 0, a_, b_, c_, d_, e_, i : = 0 , k : = 0 , m_, p_, q_, x_, y_, u : = [ ] , v : = [ ] , w_ := []) 

Prog 

If k = 0 

k := DIM(f) + 1 
p_ := FL00R(C2-n) A Cl/4)/Jm) 
m_ := HIN(b A 2 , 10- b) 

e_ := ADDOINCO, VECT0R(1, j_, DIM(f) - 1)) 

Loop 

i_ :+ 1 

Loop 

p_ := N EXT_P R IH E ( p_) 

If DACOBICn, p_) = 1 exit 
a_ := p_ A 2 

b_ := SQUARE_R00T(n, p_) 

b_ :+ (n - b_ A 2)- INVERSE JI0D(2-b_, p_) 

b_ := HO DS(b_, a_) 

c_ := (b_ A 2 - n)/a_ 

q_ := INVERSE_MOD(p_, n) 

x_ := -m 

Loop 

y_ := a_-x_ A 2 + 2-b_-x_ + c_ 

d_ := fp 

Loop 

d_ := GCD(y_, d_) 

If d_ = 1 exit 
y_ :/ d_ 

#4: If ABS(y_) < m_ 

Prog 

(This is not the complete code of the program getrelQ. Please refer to the DERIVE file, Josef) 


User=Si mp(User) : 
#5: 

User=Simp(User) : 


getrel () = true 


22789 


10 2 1 


23065 0 4 3 0 


[x_rel , E_rel ] = 23052 , 1002 


23039 


14 2 0 


User=Simp(User) : 


23026 J L 0 0 1 2 


0 2 1 


SOLVE (MOD (2 2789 , n) = M0D(- 2 -3 -5 , n)) = true 


User=Si mp(User) : 


User=Si mp(User) : 


2 4 3 0 

SOLVE (MOD (23065 , n) = MOD(2 -3 -5 , n)) = true 


2 0 0 2 

SOLVE (MOD (2305 2 , n) = MOD(- 2 -3 -5 , n)) = true 


User=Si mp(User) : 


User=Si mp(User) : 


2 4 2 0 

SOLVE (MOD (23039 , n) = M0D(- 2 -3 -5 , n)) = true 


2 0 12 

SOLVE (MOD (23026 , n) = MOD(2 -3 -5 , n)) = true 
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The next point on our agenda is to find subset's of all the vj , such that their product 
is a perfect square, i.e. has a vector representation with only even numbers, which is the 
null vector 0 of suitable size when considered mod 2. In other words, we must find vec- 
tors x such that x*E_rel = 0 mod 2. Those vectors form a subspace of the correspond- 
ing vector space, which can be given by a basis. 


User : 

kernel (c, e_, i_, j_ := 0, p_, t_, u_, v_, w_) := 

Prog 

c := MOD(c , 2) 
t_ := DIM(c’) 

e_ := IDENTITY_MATRIX(DIM(c)) 
u_ := {1, . . . , DIM(c) } 

Loop 

Loop 

j- :+ 1 
If j_ > t_ 

Prog 

d := SELECT(Z(c|k_) =0, k_, 1, DIM(c)) 
RETURN VECTOR(e_j,k_, k_, d) 

#12: v_ := SELECT(cik_|j_ =1, k_, u_) 

If v_ t {} exit 
p_ := FIRST (v_) 
u_ := u_ \ {p_} 
v_ := v_ \ {p_} 

Loop 

If v_ = {} exit 
i_ := FIRST (v_) 
w_ := M0D(cj.i_ + cj,p_, 2) 
c t i _ := w_ 

w_ := M0D(e_j,i_ + e_j,p_, 2) 
e_j.i_ := w_ 
v := REST (v_) 


User=Si mp(User) : 
#13: 


kernel (E_rel ) 


0 0 110 
0 10 0 1 


User=Si mp(User) : 


#14: 


0 0 110 
0 10 0 1 


10 2 1 
0 4 3 0 

10 0 2 
14 2 0 

0 0 12 


2 4 2 2 

0 4 4 2 


The computation above shows that we are very close to our ultimate goal of finding solu- 
tions of x / '2=y / '2 mod n, by simply combining the components of x_rel and the row vec- 
tors of E_rel in the way indicated above by the l’s in the vectors of the kernel k, as 
well as of the the matrix (k E_rel)/2. See the examples below for details. 
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User=Si mp(User) : 


#15: SOLVE MOD M0D(x_rel -x_rel , 

L L 3 4 


2 J 

n) , n 


2 112 A 

= MOD((- 2 >3 *5 ) , n) = true 


Si mp(User ' ) : 


2 2 

SOLVE (MOD (24792 , n) = MOD((-60) , n)) = true 


User=Si mp(User) : 

#17: [GCD(24792 - 60, n) , GCD(24792 + 60, n)] = [229, 109] 


User=Si mp(User) : 


#18: SOLVE 


Si mp(User ' ) : 


(mod(mod(x_ 


rel -x_rel , n) , n| = MOD((2 -3 -5 
2 5 J 


2 2 

SOLVE (MOD (24454 , n) = MOD(180 , n)) = true 


1 2 1 

) , n) = 


User=Si mp(User) : 

#20: [GCD(24454 + 180, n) , GCD(24454 - 180, n) ] = [109, 229] 


Ok, it's high time we put all this together into one single routine that will do all the 
steps above for us. Here, o is an offset as to the number of relations with regard to the 
size dim(f) of the factor base (usually the default value o=l will do) and the memory 
consuming matrices x_rel and E_rel will be purged by default after their use. (This may 
help to keep the file small, if you want to save the session after applying qs() to a huge 
number n, as global variables are also stored.) I also included in a textbox the output of 
all the divisiors of n found by this routine, trivial as well as nontrivial ones, so you might 
get a better feeling for the overall success rate of this routine. 


User : 

qs(n, o := 1, dispose := true, d_, k_, x_, y_) := 
Prog 

setup(n) 

getrel (DIM(f) + o) 
k_ := kernel (E_rel ) 

If k_ = [] 

RETURN "Too small number of relations!" 
x_ := FIRST(EXP(k_- LOG(x_rel )) ' ) 

#21: y_ := EXP(k_- E_rel /2 • LOG(f)) 

If dispose 

[x_rel :=, E_rel :=] 

d_ := VECTOR(GCD(z_, n) , z_, x_ + y_) 
DISPLAY(d_) 

d_ := SELECT(1 < z_ < n, z_, d_) 

If d_ = [] 

RETURN 1 
MIN(d_) 

[229, 109] 


User=Simp(User) : 
# 22 : 


qs(24961) = 109 
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Well, this was a rather small number, hence let's conclude by trying a bigger one, viz. 
the notorious Mersenne number n=2^67-l, which Mersenne himself thought mistakingly 
to be prime. Needless to say that for the numbers if this size the overhead of the 
quadratic sieve is still to big to be a real competitor to the other factoring methods we 
have been discussing so far. All the same, if you have a look at the huge matrix below 
you might still wonder how Derive can be so fast, namely less than 3s on my PC! As al- 
ways, have fun when trying out examples of your own! 

User=Si mp(User) : 

67 

#23: setup(2 - 1) = true 

User=Si mp(User) : 

#24: [n, DIM(n) , b, DIM(f) , m] = [147573952589676412927, 21, 312, 31, 94] 

User=Si mp(User) : 

#25 : getrel () = true 

User : 

#26: [x_rel, E_rel ] 

Si mp(#26) : 


T 85138228182631575566 iril22100110020000 

25432000691266516464 2161210110000000 

27131959020989732516 1021222010000100 

91834232634966827117 0121020200000000 

21205539602124754815 1122000110000010 

138631399311064240574 0251110000100020 

134120451326901583109 0 1 0 0 0 0 1 0 1 0 1 0 1 0 0 0 

17080547131876123141 0111111010000000 

61343671582884063679 1131002011001001 

(This is again only a part of the output. Please refer to the file, Josef) 

User : 

#28: kernel (E_rel ) 


Si mp(#28) : 

'000110100101110100111010 
#29: 001010101011101100101110 

.010101110100110100001010 


11010000 

11000100 

00000011 


[761838257287, 193707721, 147573952589676412927] 


User=Si mp(User) : 
#30: 


67 

qs(2 - 1) = 193707721 
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These are the abstracts of the special session "Applications and 
Libraries development in Derive" 

G. Aguilera, J. L. Galan, M. A. Galan, Y. Padilla, P. Rodriguez, 

University of Malaga 

Generating random Samples from continuous and discrete distributions 

In this talk we will introduce the utility file RandomDistributions.mth. This file has been 

developed for generating random values from main continuous and discrete distributions. 

The programs contained in the file can be grouped within the following blocks: 

• Random values from uniform distribution: the program RandomUniform returns 
an uniform random sample between 0 and 1. This program is the base of the other 
generations. Different algorithms have been used to develop this program in order to 
improve the Derive's build-on function to generate samples from a continue uniform 
random distribution. 

• Random values from discrete distributions: some generic algorithms for any 
discrete distribution have been implemented as well as specific algorithms for some 
discrete distributions (Uniform, Poisson, Binomial, Geometric, Negative Binomial, ...). 

• Random values from continuous distributions: specific algorithms for main con- 
tinuous distributions have been implemented (Uniform, Exponential, Normal, Log- 
normal, Cauchy, Chi-square, Student's t, F, Gamma, Beta, ...). 

• Graphical approach: a program to plot the obtained samples together with the 
density function has been developed. With these drawings we can check graphically 
if the generated samples fit the distributions. 

The use of this utility file is useful for simulating any process which follows a specific dis- 
tribution. 


Michel Beaudin, Ecole de technologie superieure, Montreal 
Another Look at a Trusted Mathematical Assistant 

From the DERIVE user manual (version 3, September 1994), we can read the following: 
"Making mathematics more exciting and enjoyable is the driving force behind the devel- 
opment of the DERIVE program". In this talk, we will try to show how some mathemati- 
cal concepts, studied by engineering students at university level - differential equations, 
multiple variable calculus, systems of non linear equations -, can be easily illustrated by 
DERIVE. Some will object that any other CAS could do the same: well, this is probably 
true but, according to us, not as quickly and naturally: "To accomplish this DERIVE not 
only has to be a tireless, powerful and knowledgeable mathematical assistant, it must be 
an easy, natural, and convenient tool". Consequently, time can be spent to prove some 
theorem or formula and the computer algebra system helps to reinforce the mathemati- 
cal concepts. Our examples will also make use of new features added in the latest ver- 
sion of DERIVE (version 6.10 released in October 2004); features that were not exploited 
as should be - DERIVE has never been enough documented. But we are still convinced 
that Derive 6 was "Far too good just for students" (http://www. scientific- 
computing. com/scwma rapr04derive6.htm I ). 
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Josef Bohm, ACDCA & Technical University of Vienna 
CAS Tools for Exercising 

There is no doubt that even in times of CAS a certain amount of manipulating skills in 
various fields of math education is still necessary. 

Students need more or less exercising for mastering expanding and factorizing expres- 
sions, finding GCD and LCM, solving triangles, applying differentiation and integration 
rules to name only some of the fields where training of skills might be useful. 

We present a respective library developed in Derive reaching from set theory to calculus, 
which can support the students (and teachers as well) offering random generated prob- 
lems together with the respective solutions. 


A. Garcia, F. Garcia, G. Rodriguez, A. de la Villa, 

Univ. de Madrid, Univ. de Salamanca, Univ. Pontificia Comillas 
Toolboxes with DERIVE 

The European Area of Higher Education implies a profound change in the Spanish univer- 
sity. We are heading towards a competency-based teaching and a learning model with 
greater autonomy for the student, who becomes the centre of the educational model. 
New technologies can play an important role in this new scenario. This paper suggests a 
"new" possibility in the use of new technologies: The design of a "toolbox" which could be 
used later on by the student when needed in other subjects of their curriculum or in their 
careers. This "toolbox" will contain instructions, which may either be suggested by the 
teacher or developed by the student as being interesting and convenient so as to explore 
the mathematical concepts associated with them. 

There are two examples of possible "toolboxes." In the first one, instructions are related 
to geometrical aspects of the plane which are studied in high school. The instructions 
are concerning with drawings, equations of different geometric objects and some dis- 
tance between the quoted objects. 

In the second one, a toolbox is built for the subject Calculus of a single variable in 
engineering studies. In this toolbox there can be instructions about complex numbers to 
express them in their different expressions, calculate their power, roots, etc. It also in- 
corporates the instructions that enable to analyze the concepts of differential calculus 
such as the tangent, the study of increasing and convexity, etc. 

About integral calculus several tools can be included for the geometric applications. Also 
procedures for the trapeze method and the Simpson method are implemented. 

In the study of approximate methods of solving equations we can encourage students 
touse algorithms NEWTON and FIXED_POINT, which DERIVE has incorporated, and de- 
sign a procedure for the method of the bisection of the interval. 
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E. Roanes Lozano, F. A. Gonzales Redondo, 

Univ. Complutense de Madrid 

An Implementation on the Mayan numbering system in DERIVE 

The Mayan number system is a base 20, positional (to be read from top to bottom, not 
from left to right) system that makes use of a symbol representing zero. It has slightly 
different variations when used for counting days (in religious and astronomical contexts). 

Therefore, 20 symbols are needed to represent 0,1,2,..., 19. Of these, the positive ones 
were represented using dots (the value of each dot is 1), and horizontal segments (the 
value of each segment is 5), while the zero was denoted by a shell (see figure below). 

If a number is greater than 20, the symbols corresponding to units, twentieths, 400's, 
8000's... are stacked from bottom (units) to top in pure base 20, while in our decimal 
system, the different orders correspond to tenths, hundreds, thousands, etc. 

We can choose the input and output bases in DERIVE in Options > Mode Settings from 2, 
8, 10 and 16. Therefore we have implemented a procedure that allows to convert num- 
bers between any bases, and that returns the output in (row) vector style. We have im- 
plemented another procedure that builds the 20 Mayan symbols for 0,1,2,... 19 (making 
use of the DISPLAY command). Finally, another procedure (denoted Maya), that uses the 
previously mentioned procedures, converts any number from base 10 to base 20 and 
represents it in the Mayan numbering system. These procedures only make use of the 
standard DERIVE commands. 


We do not know of any other similar implementations in CASs, apart from [1] (that uses 
special facilities for inserting graphics) and the similar [2], by the same authors of this 
paper, but written in Maple. 

We believe this is an interesting example of synergy among different branches of knowl- 
edge (Mathematics, History of Mathematics and Computer Science), that can increase 
the interest of the students for different topics. 


K. Schmidt, Schmalkalden University of Applied Sciences 

Making Life in an Introductory Linear Algebra Curse Easier with DERIVE 

In teaching linear algebra we have to deal with the following problem: while the level of 
the mathematical skills which are required to work with examples is generally low (stu- 
dents only need to add, subtract and multiply), the number of calculations is usually 
large. Therefore, working with examples is time-consuming and error-prone if done by 
hand. Students get tired quickly and lose their interest in this increasingly important area 
of mathematics. The faculty therefore decided to move the introductory linear algebra 
course from the classroom to the PC lab, and acquired a special DERIVE license that al- 
lows its use on all the PCs the faculty owns, and also on the private PCs of the students. 

A utility file was then developed to facilitate teaching by providing functions for the com- 
putation of zero matrices and vectors, matrices and vectors of ones, as well as idempo- 
tent and orthogonal matrices, "just-in-time" whenever they are needed during the 
course. The utility file also contains functions that test if a given matrix is symmetric, 
idempotent or orthogonal. 
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Jim FitzSimons 

DERIVE or Maple can not factor this, x 5 - 5x + 12 
I have been told a radical root exits for this polynomial. 

Richard Schmitt [schmittrichard0yahoo.com] 

Suggestion: 

"The multi-valued nature of inverse functions" by Jeffrey 
Available at 

http : //www . apmaths . uwo . ca/~dj ef frey/offprints. html 

I have a question, too: 

Is the Lambert W function available in Derive 6? Tha nk you. 


Jim FitzSimons 
Here is all I have: 


#1: Preen si onDi gits 64 

#2: Notation Dig its := 64 

#3: InputMode := Word 

#4: fl(y) := y -t- LN(y) 

f6(x) uses iteration. 

x - y 

#3: f6(x) := ITERATE (e , y, 0,27, 10) 

f8(x) uses Newton's method. 

f y-CLN(y) - X - 1) 

#6: f8(x , yO) := ITERATE , y, yO, 32 

c y + 1 

f9(x) := 

If X = -1 + TT- i V X = -1 - 7T- L 

#7: -1 

fS(x, IF(RE(x) < -1, f6(x), 0,3)) 

#S: Wri'ghtomega(z) := f9(z) 

Test the results. 

#9: al := VECTOR ([x, f9(x)], x, -4, 4, 0,03) 

VECTOR (fl(y ) - y , y , al) 

#10: 2 1 


#11: a2 := VECTOR ([x, f9(LN(x))], x, -4, 4, 0,03) 


I must admit that I have not heard about Lambert's W function before. I found some useful 
information in CRC Concise Encyclopedia of Mathematics: 

Lambert's W-function is the inverse of the function f(W) = W e w . 14/(1) is called Omega Con- 
stant and can be considered a sort of “Golden Ratio” of exponentials since 

e -w( 1 ) 


= 14/(1). 
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I checked Jim's procedure calculating the Omega Constant: 

#12: [1, f9(LN(l))] 

#13 : [1, 0, 56 714 3 29 04 09 78 38 7 299996866 2 2103 55 549 7 5 3815 78 7186 512 5 0813 5131079 2] 

x 

#14: N SOLVE (1 = x-e , x) 

#15 : x = 0 , 56 714 3 29 04 09 78 38 7 299996866 2 2103 5556 71469 246 5 514 2 046156 5 5 3103 3 2 39 

- 0.56 714 3 29 D4 09 78 38 7 299996866 2 2103 5 5 549 7 5 3815 78 7186 512 5 0813 5131079 2 

#16 : a 

#17 : 0. 5671432904097838729999686622103555497538157871865125081351310792 

In MATHEMATICA is this function im- 
plemented as ProductLog[z]: 

ln[1]:= Product Log [ 1 ] 

0ut[1]= W{ 1) 

ln[2]:= Product Log [1.0] 

0ut[2]= 0.567143 

ln[4]:= N [Product Log [ 1 ] , 20 ] 

0ut[4]= 0.56714329040978387300 ... , , . 

The graph shows f(W) and its inverse - Lambert s 

W-function. 

Among many other papers you can download another of D. Jeffrey's papers from the above 
mentioned website: 

R.M. Corless, H. Ding, N.J. Higham, D.J. Jeffrey: The solution of S*exp(S)=A is not always the Lam- 
bert W function of A. In: "ISSAC 2007", Editor C.W.Brown, pp 116-121, ACM Press, 2007 

Nils Hahnfeld, Virgin Islands 

Hallo Josef 

How can I convince the TI-89 that the 2x2 matrix [3 2; -2 -1] does not have complex Eigenvalues but 
twice the Eigenvalue 1. Do you know any tricks? Or is it necessary to take the imaginary part 
< 0.000001 as 0? Would be nice if the next OS would perform better . . . 

Three more questions: 

- Is it possible to check the language setting of the device within a program? 

- Is it possible to present an Integral in pretty print without performing an integration? 

- When will it be possible to have a dialogue in TI-N spire programs? 

Best regards 
Nils 
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These are the matrices which were sent. The first one gives the strange Eigenvalues: 



DERIVE shows the correct Eigenvalue, but does not consider that 1 is a double Eigenvalue. 


I try to solve the characteristic equation. This seems to work similar to the DERIVE proce- 
dure: I receive the Eigenvalue 1 but not its double occurency. 



Now I check for a very small imaginary part (this is what you proposed to do). Multiple Eigen- 
values are presented. 
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iri-?m)T fet t f:-=t y fh^ y f? y y i 

| ▼ f— |R 1 geb ra Ca 1 c |0 then- |Pr gn 1 0 |C 1 ean Up | 

■ ny_eigul 

■ ny_eigvl 

0 2 -1 
2 2 5 

a 

-2 -2 2T 

4 5 0 

-4 -2 3 

000 2.000 3.000} 

£ 1 . 500 + 1 . 936 ■ i 1 . 

500- 1. 936 -i 3.00> 

...[[-2.-2.2K4.S. 


MAIN 

F;hD EKACT 

DE 12^30 


Her is another example with a double solution of the characteristic equation: 


friTmnY FT? Y F3t Y fht y fs y nr? 
f— |FHgebra|CalG|other|PrgmIo|ciean Up 



r-2 -2 2 ] 


■ eigUl 

4 4-2 



.-2 -1 3 . 

U. 000 2.000 2.000} 


1 ny_eigu 


-2 -2 2 
4 4-2 

-2 -13 


x = 1 or x = 2 


my_eigul( [ “2, “2,2 ;4„4, “2 ; “2 


DE 1E/30 


ITi-THOT T F3 t T FHt T F? Y Fb^ Y "| 

| t f— |R 1 geb ra |C a 1 c |0 ther |pr gn 1 0 |C 1 ean Up | | 


£i 

000 2.000 2.000} 


-2 -2 2 T 


■ ny_eigv 

4 4-2 

x = 1 or x = 2 


-2 -13 ... 



[T-2 -2 2 


■ ny_eigvl 

4 4-2 

£1 2 2} 


.-2-13 


...vlCl-2 

„ “2,2 ;4,4 

. “2 : “2 . - 1 .3] ) 

MAIN 

RAD EKACT 

DE 15^30 


I must admit that I don't have any idea why the Tl delivers so strange Eigenvalues for your 
2x2 matrix. (Ti-NspireCAS shows the same behaviour!) 


To your other questions: The language cannot be checked. 

It is not possible to have the Integral character within a text. 

We hope that future versions of Nspire will enable dialogues. 

Josef 


Big Air in Summertime 

Tania Koller sent a model of a jumping snowboarder (as background picture in DERIVE). This is a 
very welcome occasion to cool off on hot summer days. Many thanks to Tania and her students. 



f(x) := - 0.58 51689 596 -x + 0. 8264316376 -x + 2.784157724 



